I ^ / / 


\ 


NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 


A METHOD FOR ESTIMATING HEAT REQUIREMENTS FOR 

ICE PREVENTION ON GAS -HEATED HOLLOW 

PROPELLER BLADES 

By V. H. Gray and R, G. Campbell 

Flight Propulsion Research Laboratory 
Cleveland, Ohio 


TECHNICAL NOTE 



2 


No. 1494 



Washington 

December 1947 


NATIONAL. ADVISORY COMlvUTTEE FOR AERONAUTICS 


TE3HNICAL NOTE NO. 1494 


A f.'IETEOD FOR INTIMATING HEAT REQUIREIvlENTS FOE ICE 
PREVEfITION ON GAS -HEATED HOLLOW PROPELLER BLADES 
By V. E. Gray and R. G. Campbell 


SUt^IARY 

A detailed method is presented for determining the temperature 
and flow of heated gas necessary for ice prevention of hollow pro- 
peller blades in fllgiit and icing conditions. The propeller, blade 
is analytically divided into a n’-imber of short radial segments, 
which are successively treated as separate heat exchangers. Exprea 
sions for the total external and Internal heat transfer are com- 
bined to determine the surface temperatures of each segment. The 
thermodjTiamic steady -flow equation is given for the internal gas- 
flow process and expressions are obtained for the radial variations 
of gas temperature and pressure within the blade. For a given 
initial gas temperature in the blade shanlc cavity, the minimum gas 
flow is determined, which will provide surface temperatures of at 
least 32° F everp/here on the heated portion of the blade. 

An expression for the required heat -source input to the gas 
is included and a formula is given for calculating the required 
blade-tip nozzle area. 

A discussion is included of the Indicated benefits to be 
derived from certain alterations of the blade internal flow 
passage. . , , . 

INTRODUCTION 

Thermal prevention of- the formation of ice on aircraft pro- 
pellers has been previously investigated by means of electrical 
and hot-gas heating. The heat requirements associated with the 
use of external, electrically heated rubber blade shoes are 
analyzed in reference 1. Fli^t investigations of thermal sys- 
tems using hot gas within hollow propeller blades have been 
reported in reference 2, and by Palmatier and'Brigham of the 
Curtiss -Wright Corporation. 
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The present analysis was made at the WACA Cleveland labora- 
tory to provide a method for predicting the hot -gas flow and the 
Initial gas temperature required for satisfactory ice prevention 
on hollow propeller blades . This method involves a system of 
progressive approximations considering successive radially dis- 
posed blade segments as separate heat exchangers. Account is 
taken of the heat required to raise the temperature of the inter- 
cepted free water to the blade-surface temperature, the cooling 
effect of the evaporation of water, and the kinetic heating of the 
external -air boundary layer. The physical changes of the hot gas 
in flowing through the hollovr blade are determined in order to 
find the required tip nozzle area. The detailed analysis is 
applied to a typical propeller blade for flight at two assumed 
operating conditions at the same icing conditions to illustrate 
the step-by-step procedure and to demonstrate typical results. 

A method of modifying the blade internal passage is suggested 
whereby the heating requirements may be reduced for a given 
application of the hot -gas meth-’d of preventing ice on typical 
hollow propellers. 


SYl'fflOLS 

The following symbols are used in equations taken from 
references : 

A heat -transfer area, (sq ft) 

Aj^ cross-sectional area of blade tip nozzle, (sq ft) 

Ap cross-sectional area of blade Internal -flow passage, 

(sq ft) 

blade-section lift coefficient 
c blade chord, (ft) 

c specific heat of gas for polytropic process, 

(Btu/(lb)(°F)) 

c^ specific heat of gas at constant pressure, (Btu/(lb)(°F)) 

Dp diameter of cylinder whose radius is equa]. to leading- 

edge radius of given blade section, (ft) 

Dh hydraulic diameter, (4 A^/P), (ff) 
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f 

g 

H 

h 


J 


k 

L 

I 

M 

Mq, 

m 

Nu 

n 

P 

Pr 

P 

Pv 

% 


Qj 


change of radial kinetic energy per pound of gas in flow- 
ing throu^ radial "blade segjsent^ (ft-l"b/l"b) 

"base of natural logarithms 

flow energy dissipated "by friction per po'uid of gas in 
flowing throu^ radial "blade segment;. (ft~lh/l"b)- 

friction coefficient 

acceleration of gravity; 32,2 (ft/sec^) 

rate of heat transfer per unit area, (Btu/(hr)(sq ft)) 

convective heat -transfer coefficient; (Btu/(hr) (sq ft)(°F)) 

mechanical eq.ulvalent of heat, (ft-l"b/Btu) 

thermal conductivity; (Btu/(hr) (sq ft) (°S’/ft) ) 

latent heat of evaporation of water; (Btu/lb) 

radial length of blade segment; (ft). 

rate of interception of water, (lb/(hr)(si ft)) 

rate of evaporation of water, (lb/(hr)(sq ft)) 

liquid-water content of ambient air, (grams.'i/cu m) 

Nusselt number, (hS^A ^S.(./k) 

exponent in polytropic process, pv’^ = constant 

perimeter of blade internal -flow paasago, (ft) 

Prandtl number, (3600g) 4 Cp/k 

absolute static pressiire, (Ib/sq ft) 

pressirre of saturated water vapo.r, (Ib/sq ft) 

heat input to gas flow per blade from hea.t source, 

(Btu/hr) 

heat escape per bl.ade at tip nozale, (E"'ui/hr) 

total heat added per blade to internal gas flow from heat 
source and propeller work, (Btu/hr) 
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%-2 


a 

R 


•®® 8,2 


r 

Sjji 


s 


T 

t 

At 

U 


1X4. 


V 


V 


av 


Vr 


hea.t transferred througli blade metal from gas flowing 
throu^^i radial blade segment, (Btu/hr) 

net heat loss per pound of gas in flowing through radial 

blade segment. *^tran3 „ ^1-2 , (Btu/lb) 
w j ^ 

dynamic pressure of ^bient air relative to blade sta- 
• tion, (Ib/sq ft) 

gas constant (for air = 53.3), (ft-lb/(lb) (®R) ) 

Reynolds number baaed' on blade-section chord 

Reynolds number based on hydraulic diameter 

Reynolds number based on laminar boundary-layer thiclmess 

radivis at any propeller -blade station, (ft) 

surface total length of heated blade section in chordwise 
direction over both camber and thrust face, (ft) 

surface distance from blade-section stagnation point to 
any point on heated surface in chordwise direction, (ft) 

absolute total temperature, (^R) 

static temperature, (*^F) 

temperature rise at surface due to boundary -layer 
friction, (°F) , ’ . 

internal energy of gas, (Btu/lb) 

radial velocity of internal gas relative to propeller 
blade, (ft/sec) 

tangential velocity of propeller blade at given radius, 
(ft/sec) 

air velocity, (ft/aec) 

average air velocity over camber or thrust face o.fc. blade, 
(ft/sec) 

resultant or helical velocity of propeller blade at any 
radius, (ft/aec) 
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V specific Tolume of internal gas, (cu ft /lb) 

^prop increment of enthalpy due to worJc done by rotating pro- 
peller on pound of gas in flowing throu^ radial blade 
segment, (ft-lb/lb) 

W ^_2 net compression work done on pound of gas in flowing 

throu^- redial blade segment, ^^prop " (ft-lb/lb) 

w rate of internal gas flow per blade, (lb /hr) 

X Hardy's evaporation factor 

X distance from l.eading edge measured along chord, (ft) 

Z exponent of Pr in determining kinetic temperature rise, 

1/2 for laminar flow and l/3 for turbulent flow 

a angle of attack of blade section, (deg) 

7 ratio of specific heats 

5^ thicloiess of laminar boundary layer on blade surface, (ft) 

5^ heat -transfer length of turbulent boundary layer, (ft) 

6 momentum thickness of boundary layer on blade surface 

at transition, (ft) 

t turbulent boundary- layer parameter 

t-ti;. turbulent boundary-layer parameter at point of transition 

A propeller drive -shaft terpue increment due to internal 

gas flow, (lb-ft/(lb/sec) ) 

p absolute viscosity, ( (lb) (aec)/(sq ft)) 

I? kinematic viscosity, (sq ft/sec) 

p density, ((lb) (sec^)/ft^) 

T static temperature, (°K) 

Ty arithmetic average of static temperatures of ambient air 

and blade surface, (°E) 

angular position relative to stagnation point on blade- 
section leading-edge circular arc, (deg) 


0 
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\1/ angle of Impingeinent of water droplet on Made surface^ 

(deg) 

d) propeller rotationa] speed, (radians /sec ) 

Subscripts : 

0 ambient atmospheric conditions 

1, 2 Interna] gas at entering and leaving ends, respective I7, 
of a given radial segment 

a external air side of propeller blade 

av average 

b outer edge of extemal-air boundary’- layer 

d datum temperature (for determining beat -transfer differ- 

entials) 

f final gas conditions (exit from last radial blade segaent 

prior to tip nozzle) 

g internal gas side of propeller blade 

1 Initial gas condition (inlet to first radial blade 

segment) 

m mean value between points 1 and 2 

s external blade surface 

w occurrence of condensation or evaporation of water 

I 

Primed symbols denote conditions that are changed because of 
altered blade interiors . 


METHOD FOB DETERMIWIWG DSSIGW REQUIREI'IEI'ITS 
Description of Solution 

The determination of the req^uired gas -flow rates for a hot- 
gas type of anti -icing system on a given aircraft propeller 
involves a tedious series of calcul.ations . The information 
available at the present time is insufficient to determine 
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q,uick].y the maximum req^uirementa for a given propeller operating 
through a range of flight and icing conditions. Individual calcu- 
lations must be made for a complete series of selected, critical 
flight and icing conditions. The blade-t:.p nozzle area and the 
heat-source capacity are then calculated for the maximimi reqxiire- 
ments so determined. 

The degree of approximation obtained through the use of the 
subsequent analysis will remain luacertain until more complete and 
reliable experimental data are available on propeller -blade 
external and internal heat -transf er coefficients^ kinetic heating 
of wet air, heat losses due to evaporation of water, the extent 
of heat coiiduction in the blade meta,l, and the determina.tiOn of 
minimum standards for ice prevention. 

A schematic diagram of the heated -gas flow throu^ a typical 
hollow propeller blade is illustrated, in figure . The gas is 
first heated by a suitable heat source, enters a stationary 
transfer manifold, and then passes throu^ a collector ring and 
cuffs fastened to the propeller hub and into the blade shanks. 

The hot gas then flows radially outward throu,^ the hollow blades, 
which may have internal radial partitions (fig. l), to discharge 
nozzles in the blade tips . 

The method suggested starts with the selection of a number 
of operating conditions (propeller rotational speed, airspeed, 
density altitude, and aiabiont-air temperature) characteristic of 
taxiing, talce-off, climb, cruise, and. maximum speed with the 
determination of the corresponding liquid-water contents from 
the recommendations shown in figure 2. The required internal 
hot-gas flow is then calculated for each condition using an 
initial hot -gas temperature as high as structural considerations 
will permit . 

For some flight conditions the required internal gas flow 
obtained in the calculations using the maximum initial gas tem- 
perature will be so low that a reduction in the initial tempera- 
ture and an increase in the gas flow might seem advisable. In 
many cases, however, the maximum heat requirements will probably 
be met only by utilizing the highest possible initial hot-gas 
temperature because of limitations of the propeller in pumping 
the gas. The use of high gas temperatures and low gas flows has 
the additional advantago of keeping the pumping power at a 
minimum. 

The required gas flow, which is defined as the rate that will 
maintain the heated portion of a blade surface at temperatures at 
or above 32*^ F, must be computed by trial. For conventional 
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propeller-tlade designs^ the point of lo’irest calculated surface 
temperature will usually he at the stagnation point of a radial 
station near the huh. Under severe icing conditions, there is no 
great variation of calculated leading-edge surface temperat'ure 
with propeller radius at the inhoard stations, although heat con- 
ducted radially outward from the hlade huh will actually raise the 
surface temperatures at the inner radii. If the effect of this 
heat conduction from the hlade shank does not ertend too far along 
the hlade. at the leading edge, the use of the stagnation-point 
surface temperature at the Innermost hlade segment as the critical 
point for determining the required gas flow will permit more rapid 
progression with the trial values and seldom will a complete 
analysis of all segments he necessary. 


Steps in Solution 

For each set of fli^it and icing conditions the following 
step-hy-step method of solution is suggested. 


I. Divide the hlade into a convenient number of radial seg- 
ments each I long. (In the numerica], example given 
subsequently, radial segments 1 ft long were employed; 
experience may indicate a need for shorter segments of 
perhaps 6 in; in length.) Tabulate the cross-sectional 
flow areas Ap, the external and internal heat-transf er 

areas Ag^ (= Sq. l) and Ag (= PZ), respectively, the 
diameters of equivalent leading-edge cylinders Dq at 
the centers of the segments, the internal perimeters P, 
and the chord longbhs c. The cross-sectional flow 
areas should he tabulated for the inlet, the center, and 
the outlet of each segment. 

II. Determine the external factors that are unaffected by 

the internal heat flow as follows : For the assumed 

propeller-blade speed, find the lift coefficient Cj 
and angle of attack a at the center of each radial 
segment by the method of reference 3 or other accepted 
procedure. With these data calculate the chordwise 
distribution of local air velocity 

A. Determine the local external air-fl3m heat-transfer 
coefficient hjj_ for the laminar regime (from refer- 
ence 4) from' 




Nu k 


( 1 ) 
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where the value of the laminar houndarj- layer thlck- 
nes3 is found from r^ferenoe 5 as 

9.17 


For the turbulent regime the heat - tranaf er coeffi- 
cient has bean given in reference 6 as 



•s/c 


a.±i 


(II) 


<l) 


( 2 ) 


■•ha = ^iy~^ (: 

^ St 

The value of the turbulent boundary -laj'-er heat- 
transfer length 5^ is slightly involved. The 
turbulent boundary r- layer parameter I, as derived 
in x-eference 7, is related to 5.j^ as follows: 




(4) 


The value of the ' parameter . denoted at the 

transition point (appendix A) , is given in refer- 
ence 7 as 


= 2.56 logg 4.075 


Vb e 


V 


(5) 


vnere 


C » . 0 # 289 &2 


(6) 


the value of being taken at the transition 

point. Succeeding values of t may be found from 
the relation 


^ ,' 6.13 ^ \ ... 

dx dx u 


(7) 


where the function f(t) is 


fit) = 10.411 


.-2 -0.3914? 

.( e 


( 8 ) 


10 


MGA TR No. 1494 


An alternate and simpler method than ^that above for 
determining hg^ (explained in appendix A) is given 

in reference 8 as 


, .O.SOYVov Pf) 

ha = 0.0562 (Ty) (^_av_iL_ j 


0.50 


for the laminar regime^ and 


= 0.534 


for the tiorbulent regime^ where 


0.80 


(S) 


( 10 ) 


V, 


av 


'K 


'1 


V 4 cos CLJ 


( 11 ) 


The + is used for the camber surface and the - is 
used for the thrust face. Inasmuch as the values of 
hg^ in equations (1), (3)^ (9), and (10) cannot be 

determined at the stagnation point, determine values 
of hg near the leading edge (appendix A) from the 
eq^uation of reference 8 


....... <v“^“ (.. I^f) ,« 

B. Determine the rate of water interception M (appen- 
dix B) from 


C. Determine the local .heat -transfer datum air tempera- 
ture tg ^ (appendix B) from 

for saturated air flowing over a wet surface. If 
the air flowing over the surface is unsjiturated, the 
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local air temperature at the outer edge of the 
"boundary laj-er ■becomes t-j^, and if the surface is 
dry the air-temperature rise in the boundary layer 
due to friction becomes 4tg_. The values of 

and t-jj are determined from the adiabatic compres- 
sion lines shown in figure 3, plotted from data of 
refex’ence 9. Fi-om the Intel’s action point of the 
ambient -air temperature tg and pressure Pq^ 

determine the proper adiabatic line,; then the local 
temperature at the outer edge of the boimdary layer 
t-jj at any point over the airfoil is found by fol- 
lowing this adiabatic ].ine parallel to these shown 
to the local pressure given by 


P- = p,. + CL 1 - 


n, 


b\ 




■; 


(15) 


The values of Atj^ and Atg^ have been developed 
in reference 10 as follows : 


At„ = 


Pr^ 


2gJc 


P 


(16) 


and 


At 







(17) 


where the saturated vapor pressures ^ 
correspond to the temperatures t. , and 

gL ^ CL 

respectively, necessitating a solution of 
ti-ial at each point. 


and 


III. In order to determine the changes in the internal gas 

flow through the blade, select the desired inlet-gas tem- 
perature which equals 'tg, i "tbe first radial 

segment, and a trial value of hot-gas flow w. 


A. Assume the temperature of the gas leaving the first 
segment t_ , . 
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1. Betenalne the internal gas-fiiin heat-transfer 
coefficient hg (appendix C) from reference 11 

, xO.3 0.8 „0.2 

,h^ = 4.1 X 10 (18) 




using 




- 


+ T 


-&X.' 


(19) 


2. Determine the mean heat -transfer datum gas tem- 
pei’ature for the segment ‘tg^d from 




(u ) Pr 


,0.33 


2gJc 


P 


(ao) 


where 

# 


%’,m 


w R g,m 


- Aj, 3600 


m 


( 21 ) 


For the first estimate of t 


g,2^ 


lot 


■g..m 


Pr- 1 J suhseguent trials let 


V, 


V ■,+ V r. 

g.l g,'^ 


g,m 


of the preceding trial. (Details for calculat- 
ing p_ p will he presented suhseguently. ) As 

the solution is approached this error will 
become very small. 


3. Determine the chordwlse distribution of surface 
temperat'ore tg. (appendix C) at the center of 
the segment from 



7 T^\l 


V*a,d + M \*0 ■> Aa 


ligAg ■+ Jv + ‘MrAl 


( 22 ) 
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whore (see appendix B) 


X = 


1 + 


Py,B ~ Py,a 

^8 " "^a,d 


/0.623 l\ 

V"Ph“7 


(23) 


and g and ^ correspond to the tem- 
peratures tg and tg^ respectiyely, for a 
wet surface. Because X varies with tg^ the 
'values of tg at each point must be computed 
■by trial. 


4. In order to determine the heat transferred 
thro’ugh the blado metal of the segaent, plot 
(tg - tg) against s/c for the chordwise 

extent of the heated surface and obtain the 
average value by graphical integration or by 
Simpson's rule. The heat transferred throu^ 
the blade metal is then 

^rans ^^g,d " ^s)av ^g'^g 

5. Determine the energy change due to friction in 
the segment F;j ^_2 per pound of gas flow (from 
reference 12 ) 


F 


m _ J_ ^ 


1-2 " 2g ' Dj^ 


(25) 


where 


f = 0.005B + i (BeD^ia)'°‘^^ (26) 


and 


Be. 


w 


" A^ l3^ g) 


(27) 


6. Determine the polytropic specific heat of 
gas c^ from 

I c, - _ *^rans "^1-2 „ j. \ /or,\ 

ftl_2 - -g= = (^g,l ■ %,2) (^8) 
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and the polytropic exponent n from 




(29) 


7. Determine the pressure of the gaa leaving the 
segment Pg^ g f 


%>,2 " ^g,l 



(30) 


where Pg 2 . estimated for a particular hot- 

gas -flow system from published experimental 
pressure-loss data on ducts and manifolds. 

Find the change in radial Icinetic energy Sp^g 

per pound of gas flow throu^ the segment from 



8. Deteriaj.ne the net compression work 

(appendix D) done on each pound of gas in flow- 
ing through the segment from 


u+ 2 - U.(^ 

>^1-2 = - Fp.2 = - F 


Fi-2 


(32) 


9. Finally, check the accuracy of the assumed 
(appendix D) from 


g,l - 'tg,2 - j ( y_i )■ nJcn 


.3 = 


7-n 


(Wp-2 - Ep_2) (33) 


and reestimate ‘fcg^2 the value obtained 

from equation (33) agrees with the initial 
estimate in step III -A. 
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B. For the subBequeut redial segments, repeat step III -A 

using the values of segment 

as tg^]_ and Pg^i; respectively, of the nerb 

segment. 

C. If tg at any point (found in step III -A -3 for the 
correct tg g) proves to be less than 32'^ F, it can 

be increased by increasing the hot -gas flow v and 
repeating step III. Usii.ally tg will be a mininium 

at the stagnation point and at a segnent near the 
hub. 


TV. From the final conditions of the gas leaving the last 
segment, determine the tip nozzle area required to 

maintain the gas flow (for unity discharge coefficient) 
from 




n 


(34) 


■0 


\i 



7-1 

Lgjc ^ 

M ' f 

J M Pq / 

[V I’D J V 




- 1 


where 


Tg,!- = * 


u. 


r,f 

2gJc 


(35) 


V. 


Calculate the follovring: The required input from the 

heat -source Qj, excluding losses in the induction sys- 


tem prior to the initial point, from 

r 

Ql = w Cp (Tg^i - Tq) - j 


(36) 


V. 


the total heat added to the gas flow Qj. from 


Qj = Qj 


• J 2g 


the heat escape at the tip nozzle from 


( 37 ) 
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(38) 


and tlie heat transferred throu^i the entire heated sur- 
face of the hlade liOt^ans 




rans 


«n - » <=p ('fg,l-Tg,f) + I 






(39) 


Thus, the effectiveness of the propeller blade as a heat 

^^trans 

exchanger is given ty — . 

tjji 


TTOMERICii 3XA14PLSS 
Flight and Icing Conditions 

The procedure for determining the design requirements for a 
gas -heated propeller has been applied to a theoretical hollow steel 
propeller having blades with central radial ribs, as shown in 
fig'ire 1. The propeller -blade -form characteristics are given in 
figure 4 and ocher pertinent data in table I. An assumption has 
been made that heated air enters the blade shanks from a colD.ector 
ring through suitable orifices and passes radially outward through 
only the leading-edge cavity in the blades, as illustrated in 
figure 1. The blade was arbitrarily divided for convenience into 
four radial segments each 1 foot long starting at the 18 -inch 
station. The radial stations for which the surface-temperature 
calculations were made were therefore at 24, 36, 43, and 60 inches. 
(See fig. 4.) For more accurate results calculations nearer the 
hub may be necessary. A greater number of shorter radial segments 
may also be needed in regions where either the blade airfoil sec- 
tion or the internal flow area changes rapidly with radius. 

For purposes of comparison, the calculations have been made 
for the two following operating conditions through the same icing 
condition: 


Condition 

A 

B 

Pressure altitude, ft 

18,000 

18,000 

Ambient-air temperature, °F 

0 

0 

True flight speed, mph 

400 

325 

Propeller speed, rpm 

1430 

800 

Liquid-water content, g/cu m 

0.4 

0.4 
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The value of liquid -water content chosen was based on the 
recoinniendations shown in figure 2. 

For convenience in determining the initial internal gas condi- 
tions at the blade shanks, an arbitrary assumption was made that 
the internal static pressure in the blade shanks would be the 
ambient atmospheric pressure plus three-fourths of the fli^it 
dynamic pressure. The assumed hot -gas flow was varied at an 
initial temperature of 500° F until surface tempei-atures of 32° F 
or more were obtained over the heated surface of the blade. 


Results 

The general results are summarized as follows : 


Condition ] 

1 A 

■ B 

Gas flow per blade, Ib/hr | 

450 

750 

Initial gas temperature, °F I 

1 .500 

500 

Final gas temperature, °F 

.349.6 

320.8 

Heat-source input required per blade, Btu/hr 

50, 550 

86,680 

Required blade-tip nozzle area, sq ft j 

0.00466 

0.0132 


The results of the calculations are given in more detail in 
table II and figures 5 to 10. For the four radial stations and 
for both simulated fli^xt conditions, the chordwlee variation of 
local velocity ratio is shovm in figure 5; the external 

heat-transfer coefficient h^ is shown in figure 6; the- local 
water interception M and evaporation rates are shown in 

figure 7; the local heat -transfer datum temperature tg^^^ is 
shown in figure 8; and the external-surface temperature tg is 
shown in figure 9. The radial distributions of internal gas tem- 
perature tg, internal gas pressure p^, internal heat -transfer 


coefficient 




stagnation point 
and B. 


S’ 

and the .external -surface tempei-ature at the 
tg are shown in figure 10 for conditions A 


The calculations of the external heat -transfer coefficients 
v;ere made using equations (1), (3), and (12). Attempts were made 
to locate accurately the external boimdary-layex- tz’ansition 
points using the methods of references 13 and 1C, and the theoret- 
ical transition points were found to be located aft of the heated 
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areas. In four cases the velocity distributions (fig. 5) have 
fairly distinct peaks near the leading edge and an a.rbitrary 
assumption was made that the surface rou^ness due to the imping- 
ing water would cause transition at these minimum-pressure points. 
This assumption appl.ies to the camber face at the 48 -inch station 
for condition A and to the thrust face at the 60-inch station and 
the camber face at the 36- and 48-inch stations for condition B. 
Ihe resulting chordwise distributions of hji, sho^-ai in figure 6 

for these cases, are similar to tile type sho'tm on the left aide of 
figure 11 and discussed in appendix A, inasmuch as in each case 
the transition point occurs very close to the leading edge. 

The surface temperatures are obtained fi’om eq_uation (22), in 
which the external heat -transfer area Aa is arbitrai’ily assumed 
to extend aft of the blade-rib center line for a distance on each 
face equal to half the maximum blade thickness . This assianirtion 
is to account approximately for the heat ti-anafer to the trailing 
half of the blade, which in reality diminishes all the way to the 
trailing edge. An unpartitioned blade would not require such an 
approximation. 


Discussiorj • 

Thin-Skin Approximation 

The largest source of error in this analysis is thought to be 
the thin-skin apnroximation used in the solution df extexnal sui'- 
face temperatures. This approximation makes use of two related 
assumpuions. First, the value of Eg depends on the temperature 
differential ('fcg,d " '*^s)^ which the temperature gradient 
through the blade metal is assumed to be zero. This ^adient and 
the error involved is usually small, except at the leadj.ng-edge 
region, for any propeller blade made of steel or metal with equal 
or higher thermal conductivity. Second, the heat balance employed 
in equation (22) is based on the assumption that the total heat 
transfer from the internal gas to the 'external air can be con- 
sidered to be distributed uniformly and to pass only noiiaally 
throvi^ suitable internal and external heat -transfer areas; that 
is, the blade metal is assumed to be infinitely thin and surface- 
wise heat conduction is neglected. The ratio of the two heat- 
transfer areas is then used to determine the ratio of the local 
internal and external heat -transfer rates. The error involved in 
this assumption is thouglit to be large inasmuch as the thermal 
conductivity of the blade metal is usually very hi^^i compared 
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vith the conductivity through the gas and air aui-face f'i.ljns. Sur~ 
facewiae heat conduction within tlie "blade metal is certain to 
reduce the temperature gradients indicated hy the variations of 
surface-film heat-transfer coefficients and changes in metal 
thicloiess. 

The two foregoing errors can he eliminated from the analysis 
if the blade segments are studied using the "relaxation" method 
developed in reference 15. Wo attempt has been made to use this 
method because of the prodigious amount of t.lmo and effort involved. 
The method and a sample of the results to be expected from it are 
given in appendix E. 


Principal Factors Affecting Heat Requirements 

From a consideration of the foregoing analysis and the numer- 
ical examples, the required internal heat flow through, hollow pro- 
pellers depends on the following factors external to the blade 
surface, which appear to have the largest effects in determining 
the local heat transfer through the metal necessary to prevent Ice. 

Ambient tempera ture . - With other factors remaining constant, 
the heat transfer through the blade surface required for ice pre- 
vention as defined herein varies directly with the temperature 
difference between the freezing temperature and the datum tempera- 
ture (ambient temperature plus kinetic increment). Even thou^ 
the assumed liquid-water content of the ambient air deci’eases with 
decreasing temperatures, the heat required to maintain a minimum 
surface temuerature of 52*^ P steadily increases as ambient tempera- 
ture decreases j but, computed on the basis of immediate evapora- 
tion of all the water that strikes the surface, the heat required 
will decrease with decreasing temperatures because the liquid- 
water content decreases. The point at which the heat requirements 
are equal when computed on both bases of heating represents the 
maximtmi possible value of heat required. This maximum value rep- 
resents excessive heating for the conditions used in the examples 
presented. 

Water impingement. - An increase in the rate of water inter- 
ception requires a direct increase in the amount of .heat required 
to raise its temperature above 32° F. A more important factor, 
when external heat requirements are considered, is that an increase 
in the rate of water interception increases the •■amount of runback, 
which increases the area of wetted surface apd "thus requires. mors 
heat to offset the evaporative heat loss. The difference in tem- 
peratures between a wet and a dry surface resulting from the same 
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internal heat flow can he ohserved' in figure 9(h), where the 
surface hecomes dry on the thrust face of hoth the 48- and 60-inch 
stations. Tlie atirface-temperature rise ahovo ambient temperature 
more than doubles its value in the two instances after the surface 
hecomes dry. 

Kinetic heating . - Kinetic heating always reduces the heat 
requirements for ice prevention because it increases the heat- 
transfer datum temperature above the ambient tempei’ature . Its 
temperature rise increases as the square of the resultant velocity 
and has a lower value in wet aii’ than in dry air. This increment 
of temperature becomes important as a saving in required heat flow 
at high resultant velocities of the order shown in the two examples 
presented. A fact of interest is that heat requirements are 
reduced according to the square of the resultant velocity (because 
of kinetic heating) and are also increased according to the square 
root of the resultant velocity (because of the heat-transfer 
coefficient). The net result is that as velocity is increased 
the heet required to maintain a sui’face temperature of 32° F first 
increases and then decreases ; the location of the maximvim heat • 
requirement depends on the conditions involved. This maximum heat 
requirement probably occurs below a velocity of 400 feet per sec- 
ond for most conditions. Thus, except for other variations (such 
as changes in propeller -blade sections), the heat requii-ements for 
propellers on high-speed airplanes are usually greatest at the , 
inboard stations. For the same reason, the heat requirements are 
reduced by increasing the rotational speed of the propeller. 

Location of transition. - The determination of heat transfer 
from aT prope5.1er blade*~de^nds to a large extent on the exact 
location of transition. As shown in figure 11, considerably more 
heat transfer results over the blade faces for which turbulent 
flow is assumed (as shown by the higher heat -transfer coefficients) 
than would be the case if laminar flow were assumed. The same 
result is shown in figure 9, in which the low surface temperatiAres 
of the four faces under assumed turbulent flow indicate a higher:’ 
rate of heat transfer. Of course, transition might be caused by 
water-film rou^ness on many of the sixrfaces for which laminar 
flow was assumed, because the velocity gradients are quite low 
over all the stations, as shown in figure 5. 

Alt itude. - The- altitude -pres sure term frequently appears in 
the analysis and Its total effect on heat requii’ements is difficult 
to evaluate. The two principal effects of altitude (other factors 
remaining constant) on the required blade-metal heat transfer are 
largely compensatory.- -With increasing altitude the heat transfer 
due to the evaporation processes increases whereas the heat transfer 
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due to external convection decreases. Increasing altitudes have a 
serious effect on restricting the internal gas flow. In order to 
offset the reduction in gas density with altitude, the gas flow w 
must he maintained hy increasing the internal gas velocity, with 
resulting increases in flow friction and pressure drop. These 
increases, coupled with the loss of ram pressure at altitizde, indi- 
cate that the required internal gas flow steadily becomes more 
difficult to maintain with Increasing altitude and that a booster 
pump may be necessary at high altitudes to attain the required gas 
flow. A pressvire boost can also be obtained by increasing the pro- 
peller rotational speed whenever possible. 


Comparison of Results of Numerical Examples 

Flight condition B required a larger internal heat flow than 
condition A for two reasons: the kinetic heating was leas, owing 

to a lower airspeed and propeller speed, and the surface area 
under assumed t'orbulent flow was greater. 

The required tip nozzle area for condition B was almost as 
large as the internal flow area at the exit to station 60, which 
indicated that the gas pressure at the tip had dropped nearly to 
its lowest limit and that a hi^ier rate of flow was almost 
impossible without augmenting the initial gas pressure or increas- 
ing the propeller speed. If this tip nozzle area is used, the 
flow will become excessive under conditions of hl^ier airspeed 
and propeller speed or lower altitude. TTie hot -gas flow could, 
however, be regulated by a throttle in the intake system. 

The required heat-source capacities for these two conditions 
appeared to be very high compared with external electrically 
heated blade shoes based on the meager experimental data avail- 
able. The high required heat -source input is, however, explained 
by the fact that internal gas -heating systems normally have no 
selective control of heat transfer and as a result much of the 
surface area is considerably overheated and further wastage 
occurs at the tip nozzle. In the examples that used this method 
of blade heating, only a narrow band at the leading-edge surface 
had temperatures near 32° F. Consequently, if a narrow atrip of 
ice 0.3 inch wide, for example, were to be permitted along the 
leading edge, the internal heat -flow requirement for station 24, 
condition B, (most critical) would be nearly halved. This 
reduction would result because the surface temperature at the 
stagnation point could then be allowed to drop approximately 
9° below 32° F and the differential (tg - tg , 3 ^) would be 
approximately one-half its original value (17!s° F), Because 
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the stagnation point is the controlling factor for the segment^ 
the inteznal heat flow could then he reduced by nearly half and 
the blade-surface temperatures would remain above 32*^ F except for 
the narrow leading-edge band. In practice^ the internal heat flow 
required therefore depends lai'gely on whether narrow strips of ice 
along the leading edges are tolerated. 


A method of reducing the required heat flow in a typical hol- 
low propeller blade is suggested in which the internal flow passage 
is altered to attain a more economical distribution of heat transfer. 
When it is assiJmed that equation (18) and the hydraulic diameter 
concept can be applied to internal passages whose perimeters con- 
tain convolutions, designs of flow passages can be developed whereby 
large reductions in heat-source capacity and pijmping -power output 
are apparent over the typical hollow-blade interior, as shown in 
figure 1. G?wo suggested alterations are shoTO in figure 12 and 
compared with the original blade section. Using primes to denote 
the altered designs, the following mathematical comparisons can 
be shown: 

For both the original and altered designs 


from equation (18). For the same inlet gas temperatures, assuming 
P = Ag and taking the ratio of a modified design to the original 

design 


Suggested Technique for Eeducing Heat Flo-^ra 






Setting this ratio equal to unity will permit evaluation of the 
flow ratio that will produce the same amount of blade -surf ace heat 
transfer as the original design, or 
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Modif ipa- j Modifica- 
tion 1 I tion 2 

0.149 I 0.152 

i 

i 

i 


" A^pgT36b0)' 
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^ 1-2 


Uv 


1.68 


(%) 


-1.32 


I 

! 

} 

i 

I 

.298 I .228 



In auimnarlzing these modifications, if the area of the inter- 
nal flow passage is reduced and the internal heat-transf e'r area is 
increased, the same total heat transfer to the external surface 
can he maintained while the intemal-gas flow is gi’eatly reduced, 

A reduction of the flow lowers the internal velocity and the fric- 
tion loss and thus raises the available internal pressure, \7hen 
the flow is reduced, more heat is removed per pound of gas flow 
and, consequently, the gas temperature falls more rapidly and 
reduces the heat loss at the tip nozzle. Care should be taken 
with the internal heat -transfer area that some segments do not 
aboti-act such heat as to starve subsequent segments. 


In addition to the foregoing reduction of heat flow, a better 
distribution of chordwise surface heat transfer can be accomplished 
by making the greatest increase of internal heat-tz’ansfer area at 
the leading edge by the use of fins similar to those shoTO in 
figure 12. An accurate determination of the magnitude of the heat 
transfer resulting from these modifications requires that the 
whole problem be approached along the lines of the relaxation 
method as outlined in appendix E. All the indications shown lead 
to the belief that a hollow iDropeller blade can be designed with 
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modif icationa similar to those illustrated., -which vill efficiently 
provide ice pi-otection with a heat flov of feasible proportions. 


COMMENTS AND RECOMMENDATIONS 

The method presented for estimating heat roQ^uirements for ice 
prevention on gas -heated hollow propeller blades provides inte- 
grated means for determining the design requirements for ice pre- 
vention on gas -heated hollow propeller blades operating in any 
selected condition. Several of the suggested aimplif icatlons and 
other short methods can be used where conditions do not require a 
rigorous treatment. Solutions by ti-ial are occasionally required, 
which are admittedly inconvenient, but become fairly rapid with 
increased familiarity. \ 

Substantiation of several formulas, upon which this method 
is based, by experimental measurements in actual icing conditions 
is recoma.jnded . A relaxation analysis of the whole propeller 
blade should be made to determine the effect of leading-edge 
interna], ribs and. the exact heat transfer normally, radj.ally and 
chordwise throu^ the blade aietal. 


Pligiit Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, June 27, 1947. 
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APPStIDIX A 


EXTERNAL HPAT-TRANSEER COEFFICmTS 

Airfoil co e ff loj.entB based o n bou ndary layer. - Several 
methods are’~ available for calculating the local and average heat- 
transfer coeff loiontG for the surfaces of airfoils (references -i, 
16, 17, and ll). Local sui’face values for the heat-transfer 
coefficients on airfoils (eq,uation (l)) are derived in reference 4, 
primarily on the basis of Eeynolds' analog between akin friction 
and heat transfer throu^ laminar boundary layers , and an alter- 
nate form of Reyiiolds' analogy (equation (3)) is presented for the 
turbulent boundary-layer case in reference 16. These solutions 
are detailed because they take into account the variation of local 
air velocity and boundaiy- layer thickness over the airfoil surface. 
Wo account has’ been talcen either in this ana3.yais or in the refer- 
ences, however, of the unkno'svn effect on t}ie botmdary- layer heat 
transfer of the heat addition to the boundary layer caused by a 
heated surface. 

Empirical airfoi l coefficients . - The methods of reference 8 

( equations ' ( 9 ) and (lO) ) have been added as alternate eoli’tions . 

for hj^ because they offer easier solutions by their use of 
average air velocities over the faces of the airfoils. These 
equations require a sli^t but adequate approximation for use in 
this analysis when determining the average temperature Ty, 

because the surface temperatures are unknown. However, estimates 
of average surface temperatures involving errors of leas than 

1^ percent of h^ can easily be made. 

Leading-edge cylinder coefficients. - At the stagnation point 
the value of hg^ cannot be 'determined in equations (l), (3), (9), 

and (10). The suggestion is made in reference 6 that the leading 
edge of an airfoil be regarded as an isolated cylinder in trans- 
verse flow for which heat-transfer data are available in refer- 
ence ir . Wlien these data are used and Wusselt's numbei’ at the 
stagnation point is expressed as a function of Prandtl number and 
Eeynolds number, an empirical equation for the heat-treinsf er coef- 
ficient over the forward half of the equivalent leading-edge 
cylinder (equation (12)) is developed (reference 6). Figure 11 
diagremmatically illustrates the suggested technique of fairing 
the distribution of heat -transfer coefficients over the leading- 
edge cylinder into the laminar- and turbulent -flow distributions 
for the rest of the airfoil in two tj'pical configurations. On the 


26 


NACA TN No. 1494 


■fchrust face, transition is shown to occur close to the leading 
edge and the h^ curve is faired directly from the leading-edge 

values into the turbulent. -flow curve. On the camber face, transi- 
tion is shown to occur after a length of l.aminar flow and the 
curve is arbitrarily faired into the txirbulent regime. The former 
type of distribution occurs in the numerical example, in certain 
cases on both faces of the blade. 

Location of transition . - In the foregoing considerations, 
the location of the transition point must be determined. For 
increaaing local pressures, the point of transition, which may be 
assumed to be coincident with laminar separation, may be estimated 
by the method of references 13 and 14. I^or decreasing pressures, 
the suggestion is made in’ reference 5 that the Reynolds number 
based on the laminar boundary -layer tliiclmess Reg ^ point 

of transition is between the limits of 8000 and 9500. 

The location of the transition point greatly affects the heat 
transfer from an airfoil, because the heat-ti’ansf er coefficient for 
turbulent flow is considerably larger than that for laminar flow. 
No reliable method is known for predicting the point of transition 
when the airfoil is undergoing water impingement. The presence of 
water on the surface probably causes transition to occur forward 
of the point determined ,by the foregoing methods . 
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APPENDIX B 


EEPECTS OF N^ITFD SUPFACE 


Rate of vater Inter ce p tion. - The rate of water interception 
is only slightly overestimated hy assuming that the droplets 
travel in straight-line paths intersecting the hlade surfaces. 

This condition exists when the aii'foil collection efficiency 
■becomes 100 percent. Exact collection efficiencies for airfoils 
are unavailable, although fairly complete data are available for 
cylinders. For symmetrical and low-camber airfoils the use of 
existing data for cylinders whose diameters ai’e equal or related 
to the leading-edge diameters of the airfoils can be made with a 
fair degree of auproximation. Accordingly, propeller-blade-secti^on 
leadi.ng-edge cylinders, under ordinary flight conditions, have col- 
lection efficiencies varying approximately from 80 percent near 
the hub to 98 percent near the tip. The error in assijczing 
100-percent collection efficiency for the propeller sections is 
not so great in figuring 'bhe rate of wate;c interception as in 
determining the extent of wetted surfaces, as discussed in the 
succeeding section. 

The local rate of vrater interception M is thus proportional 
to the sine of the angle '1 subtended by the relative wind 
(helical speed) vector e.nd the tangent drawn to the blade surface 
at the point being considered. The value of M is given by 
equation (13). 

Surface evaporation. - Tlie effect of evapox’ation of water 
fj.-om the surface of an airfoil in increasing the heat transfer 
was first derived by Hardy (reference 13). He has presented an 
expression foi* a local evaporation factor X, which multiplies 
the local conv’-ective heat-transfer coefficient wlierevei* the sur- 
face is wetted. This factor is given by equation (23) and the 
local rate of evaporation of water from the surface is given by 


Nliere the local rate of evaporation is larger than the local 
rate of water interception, the surface will tend to become dry. 
If the blow-off of water from the sur-face is neglected, however, 
there will be runback of unevaporated water from the leading-edge 
region when the rate of intercepted water exceeds the rate of 



(40) 
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evaporation. In this analysis the surl’ace of the propeller-hlade 
section is considered dry aft of the point s/c at which 


‘s/c 


ps/c 


M d(s/c) = 




V.''- 


where tho lower limit 0 is at the stagnation point. For a dry 
Burfacej X hecomes 1 and no surface evaporation occurs. 

H eat -transfer dat^Jira air temperature. - The heat -transfer 
datum temperatvxre tg^ ^ at any point on the blade surface is the 

s^irface temperature that would be obtained if the blade were a 
nonconducting unheated body. Tho value of tg^^ is given by 

equation (14) baaed on figure 3 and equations (16) and (17). Tlie 
nature of tg ^ for a typical propeller blade operating in either 

saturated or unsaturated air is shotm in figure 13. Tl^ere ai’e tvro 
components of tg 3 ^: the adiabatic temperatirr’e t-j^ or at 

the outer edge of the boundary layer and the temperature differ- 
ential Atg or Atg ^ in the boundary layer due to kinetic fric- 
tion. The frictional temperature rise in the boundary layer is 
slightly larger in turbulent flow than in laminar flow according 
to equations (16) and (17) (as the Prandtl number is raised to the 
one-third power in turbulent flow compared with the ono-half power 
for laminar flow) . Both the components of tg ^ are affected by 

the presence of v^ater. The frictional temperature rise in a 
boundary layer that remains saturated is less than the rise in an 
unsaturated bqundary layer, due to evaporation. Likewise, the tem- 
perature at the edge of the boundary layer for a saturated air 
stream follows the saturated air adiabatic line. 

Whether or not the air is saturated at a given point is often 
difficult to deteimine and, consequently, the determination of 
tg ^ is uncertain. In this analysis the following course has 

been adopted; For a wet surface, which exists forward of the point 
where the accumulated water Intei’cepted equals the accimnulated 
water evaporated (as defined in the prev'ious 

"^a,d “ ^b,w ^^a,w 

For a dry surface, which is assumed to begin 
point 


section) 


(14) 


ab'cuptly af't of this 
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^a,d = ■tb,w + 

Accordingly, only when the ambient air ia unaaturated ia 

■*^a,d =■ ■‘^b + ^^a 


(42) 
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APPEllDIX C 


HEAT BALANCE 

Local Internal heat transfer. - The internal gas -film heat- 
transf ercoeff icient h_ is given hy equation (18), which appears 

O 

to he the moat applicable formula for detemining this coefficient. 
The equation is based on data obtained for fully developed turbu- 
lent flow in long tubes. It is doubtful, however, whether this 
equation will give satisfactorily accurate results when applied to 
flow through hollow propellers for two reasons: (l) The nonsym- 

metrical cross sections of a hollow propeller-blade passage and 
the effect on heat transfer of changes in flow area may not be 
fully accounted for in an equation baaed on flow in straigb.t cir- 
cular tubes; and (2) no account is taken in equation (18) of local 
flow variations along a plane normal to the direction of flow 
(chordwise in a propeller blade). In a rotating propeller, the 
Coriolis acceleration produces radial velocities along the leading 
edge of the internal passage that are hi^er than those toward the 
trailing edge. The local internal heat-transfer coefficient would 
therefore probably be hi^er at the leading edge than is indicated 
by the average value given by equation (18). 

When the value of hg given by equation (18) is used, the 
local internal heat transfer for a unit area can be written as 

Hg = hg ("tg^d " ■'^s^ 

where tg ^ is the mean-flow static gas temperature plus the tem- 
perature rise due to surface friction (equation (20)). 

Local external heat transfer. - The local rate of external 
heat transfer per unit area at any point may be expressed as 


Ba = l^a X (tg - ta^d) + M Cj, (tg - t^) - (4^ 

where Cj^ is the heat capacity of liquid water, taken herein as 
unity. The first term of equation (44) on the ri^t-haud side 
gives the combined convective and evaporative heat transfer; the 
second term gives the sensible heat absorbed by the intercepted 
water during its temperature rise to the surface temperature; and 
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the third term gives the heat released by the impinging water at 
the expense of its kinetic energy. For simplicity, an approximate 
expression for Hg, can be used, wherein 

Ha - (Ha ^ <‘s ‘ *a,a> 

and the error in the ax^proximation is small except at the leading 
edge where M is a maximum. 


Heat-balance equation. - In this analysis the equilibrium of 
heat transfer to and from the blade surface at any point is 
expressed as 


.. . (46) 

.. 

When equations (43) and (44) are substituted in equation (46) 
t„ can be solved for: 




Vs, a * 


h.Xt 


a, d + 




Mlt, 


0 


2gJ/_ 


'a 


+ (V + A 


S 6 


a 


( 22 ) 


Equation (46) is based on the assumption that a thin-skin 
approximation can bo used for determining the surface temperatures. 
(See section entitled "DISCUSSION.") 
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APPENDIX D 


THERMODYNAMIC ANALYSIS OF GAS FI.OW 


Aasmnptlona . - The method for determining the temperature and 
pressure of the internal gas during its flow throu^ the hollow 
propeller blade is based on the theimodynamic steady-flow energy 
equation. The internal flow is assumed to have a uniform radial 
velocity distribution at any station. The assumption is made that 
for a short radial segment the polytropic exponent n of the 
internal flow process can be considered constant. The error 
involved will vary as the radial length Z of the segments used 
in the solution. For convenience in using the polytropic rela- 
tions for gas, the enei'gy due to flow friction is assumed to be 
heat added to the gas at the expense of wort done by the gas . 

Analysis of flow process . - The thermodynamic steady-flow 
equation may be written in terms of the energies per pound of gas 
as 


This equation states that the difference between the net compres- 
sion work done on the gas and the net heat lost from the gas in 
passing through the segment is equal to the total of the differ- 
ences between the kinetic, internal, and potential energies at the 
leaving and entering ends of the radial segment. 

For a mechanically reversible process, the sum of all the 
mechanical -energy terms in equation (47) can be expressed by 


where Ej ^_2 is defined as the radial kinetic energy term by 

equation (31). The net work for a polytropic process satisfying 
the relation pv’^ = constant is given by 



El -2 + (P2V2 - PiVq) - Wi_2 


(48) 



( 49 ) 


and the flow work can be expressed as 
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Vs - J’ri “ 



(50) 


Wlaen equations (48) and (49) are equated and equation (50) is sub- 
stitu.ted, the net compression work becomes 


¥• 


1-2 = ®l-2 + - -tg,l) 


(51) 


In a similar manner, the internal -energy term of equation (47) can 
be written as 


■^(''2 - ^ <‘«,2 - * 8 , 1 > 


When equations (50), (51), and (52) are substituted, equation (47) 
becomes 


Qi_2 = : 


n-7 R / 

U-nTrr-i) J ^ 




(53) 


However, 


Q- 


1-2 


^rans 

w 


Irl . 

J 


- cjt 


%1 ■ ^6^2 


(28) 


The gas -temperature differential throu^ the segment can now be 
expressed in terms of the energy quantities by combining 
oqxiations (51), (53), and (28) to obtain 


'6,1 ■ ^8,2 


7-n 

( 7 -I) nJc 


n 


(^1-2 ■ ®l-2^ 


(33) 


In eqviation (33) the value of is given by equation (32) , 

The Increment of enthalpy work done by the rotating 

propeller on a pov’nd of gas in flowing throuj^ a radial segment is 
given by 


W. 


prop 


^,2 




2g 


(54) 


This equation can be shown as follows : The propeller-drive-shaft 

torque increment due to the internal gas flow through a blade 
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segaent is given ly the change in moment of tangential momentum, 
which may be written in terma of torque per poiuid of gas per 
second as 

A = (55) 

s 

where is the blade tangential velocity at the radius r. The 

rate of work done on a poirnd of gas can be written as 


(joA 


8 


- ^t,l^ 


8 


(56) 


because u^ = cor. For a strai^t rotating channel, the increase in 
tangential kinetic energy between two points is obviously given by 


u 


t_5_2 


u 


t,l 


28 


for a pound of gas, so the remainder of the total work represents 
the amount of external work done on a pound of gas in flowing 
through a radial segment to produce a charge in the enthalpy of the 
gas. 
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APPENDIX E 


EFFECT OF HEAT CONDUCTION THSOUGH BLADE 
METAL AT LEADING EDGE 

A study was made using the relaxation method as described in 
reference 16 to determine the effect of heat conduction in the 
blade skin on the temperatxire distribution around the leading edge. 
The original assumption that heat is transferred only in a direc- 
tion normal to the blade skin, which was assumed to be infinitely 
thin, was expected to be appreciab].y upset by the variations of 
blade-metal thickness and of the effective external heat transfer 
aroiznd the leading edge. 


Laplace's equation, which requires the reasonable assumption 
that the blade metal is a homogeneous, isotropic solid, was used 
with conditions at the Internal boundary defined by 



(57) 


and at the external boundary by 

d.0 

k = (hpX + M)(0 - 0) (58) 

where 


k 

^ • 

^x 

y 

(haX + M) 


thermal conductivity of blade metal, 
Btu/(hr)(sq ft)(^/ft) 

(tg,a - 

heat "transfer datum gas temperature, °F 

metal temperature (surface or internal), ‘^F 

distance normal to boundeiry, ft 

internal gas heat -transfer, coefficient, 
Btu/(hr)(3q ft)(^) 

effective external heat-transfer coefficient, 
Btu/(hr)(sq ft)(°F) 
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(tg,d - 

o 

heat -transfer datum air temperature, F 


The scale used was 1 inch = 0.05 inch and for convenience the 
solution was arbitrarily originated 1 inch in chord length from 
the leading edge. The metal temperatures at the starting points 
on both faces of the blade were assumed to be slightly less than 
the values of surface temperatures that were obtained in the numer- 
ical example. The inclinations of the isotherms to the surface at 
the starting points were estimated from the anticipated nature of 
the heat transfer. The values of metal temperature obtained are 
unreliable, but the isotherm and surface -temperature patterns 
obtained are believed indicative of the trends to be expected if 
the solution were to be extended over the entire blade cross 
section. 


The calculated results obtained for simulated flight and 
icing condition B at the 24-inch station are shown in fi glare 14. 
Because heat is transferred everj'Where within the metal normal 
to the isotherms, a considerable quantity of heat flows through 
the metal of the camber and thrust faces into the ma-ss of metal 
at the leading edge, as shown in figure 14. The surface tem- 
perature distribution is changed and thus the surfacewise tem- 
perature gradients at the leading edge are greatly reduced. The 
surface temperature rise at the leading edge due to conduction in 
the metal is approximately as indicated by a comparison of the 
two curves at the bottom of figure 14. 
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TABLE I - PPOPELLER DATA 


Station 

24 

7 

j 36 

48 

1 

60 

Blade radius^ r, ft 

2 

i — 

3 

4 

5 

Blade -segment radial 
lengtli^ 1, ft 

1 

1 

1 

1 

Blade chord, c, ft 

1.05 

1.06 

1.07 

0.986 

Equivalent leading-edge 
cylinder diameter, 

D^, ft 

0.027 

0.0095 

0.0057 

» 

0.0040 

External heat "transfer 
area, Aq_, sq ft 

1.26 

l.]8 

i 

! 1.15 

j 

t 

! 

j 1.06 

Internal heat -transfer 
area, Ag, sq ft 

1.13 

1.05 

1.01 

i 

0.942 

Internal hlade-passage 
perimeter, P, ft 

1.13 

i 

1 1 OS * 

1 1 

1 

1.01 

0.942 

Internal blade -pass ago 
cross-sectional area, 
A^, sq ft 

1 

0.0558 

1 

0.0315 

0.0239 

0.0173 

Internal blade passage 
cross-sectional area at 
point 1, A^^i, sq ft 

0 . 0875 

0.0375 

0.0275 

i 

0.0205| 

Internal blade passage 
cross-sectional area at 
point 2, Aq 2 ^ sg. 

0.0375 

0.0275 

0.0205 

0.0141 

- 


i 

1 
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TABLE II - CALCULATED DATA FOE NUMERICAL EXAMPLES 



Condition A 


Condi 

1 

tion B ■■ 

Station 

24 

36 

48 

60 

24 

36 

48 j 60 

^ 1 — 

Blade-section lift coef- 
ficient, C 2 

0.413 

0.550 

0.583 

0.474 

0.570 

0.597 

0.553 

0.436 

Blade-section angle of 
attack-, a, deg 

2.6 

2.8 

2.9 

2.3 

5.0 

3.2 

2.4 

1.7 

Blade resultant velocity, 
V-p, ft/sec 

658 

740 

839 

951 

505 

539 

582 

622 

Mean internal radial gas 
velocity, f^/sec 

1 

88.0 

1 

145.8 

175.5 

217.5 

159.6 

274.5 

345.9 

465.5 

Friction flow energy, 
Fi_ 2, ft-lb/lt 

15.7 

69.4 ' 

1 

I 124.8 

239.9 

46.0 

216.4 

425.9 

966.3 

Blade -metal heat transfer, 
®^rana^ 

3726 

5293 

6070 

6280 

5404 

7920 

8540 

8674 

PolytroDic specific heat, 
c^, Btu/(1^)(°F) 

0.2941 

0.3028 

0.3157 

0.3285 

0.2275 

0.2351 

0.2264 

0.1839 

Polytropic exoonent^ n 

1 

0.441 

1 

0.478 

0.525 

J 

0.563 

i 

-0.223 

-0.077 

-0.247 ! -4.51 

i ! 
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Station j 

CJ 

36 

1 

48 ! 

eo 

24 

36 

48 . 

60 

■} 

Propeller enthalpy incre- 
ment, ^^-oron^ ft-lh/lh 

1400 

2100 

2800 

3500 

438 

656 

874 

1092 

Radial kinetic -energy 
change of internal gas, 
E^_2, ft-lh/lh 

202.0 

147.9 

188.0 

401.3 

698 

608 

S55 

2575 

Tip-nozzle heat escape, 
Qq, Btu/hr 

35,250 

59,200 

1 

Summation of hlade-metal 
heat transfer, ZQ-tyans, 
Btu/hr 

21,369 

30, 538 

Effectiveness of hlade as 
heat exchanger 

0.379 

0.341 , 
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Figure I. - Schematic flow diagram for gas-heated hollow propeller 

blade. 
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Figure 2. - Variation of liquid-water content with temperature. 
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Fig. 3 
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Figure 4. - Propel I e r- b I ade- fo rm character i st i cs. 
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Fig. 5 
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(a) Condition A. 



(b) Condition B, 

Chordwise variation of local velocity ratio for two 
d i t ions. 


Figure 5 
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External heat-transfer coefficient, ha, Btu/(hr)(sq ft)(°P) 


Fig. 6 a, b 
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(a) Station 24. 
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(b) Station 36. 

Chordwise variation of externa) heat-transfer coefficient 
for two conditions. 


Figure 6 


External heat-transfer coefficient, h^, Btu/(hr)(sq ft)(^) 
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Fig. 6 c, d 
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( c) Station 48, 
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(d) Station 60. 


Figure 6 


Concluded. Chordwise variation of external heat-t ransfer 
coefficient for two conditions. 




Rate of water Interception M and water evaporation Mev» lb/(hr)(sq ft) 


Fig. 7 a, b 
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Figure 7, - Chordwise variation of rates of water interception and water 

evaporation for two conditions. 
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Fig. 
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(d) Station 60. 

Figure 7. - Concluded. Chordwise variation of rates of water inte 
ception and water evaporation for two conditions. 
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Fig. 8 
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(b) Condition B. 

Figure 8. - Chordwise variation of heat-transfer datum air temperature 

for two conditions. 
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Fig. 9 




(b) Condition B, 

Figure 9. - Chordwise variation of externaJ-sur face temperature for two 

cond it ions. 
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( a ) Condition A. 

Figure 10. - Radial variation of internal gas temperature, internal gas press- 
ure, internal heat-t ransfer coefficient, and externa I -su rface temperature at 
stagnation point. 
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(b) Condition B. 

Figure 10. - Concluded. Radial variation of internal gas temperature, in- 
ternal gas pressure, internal heat-t ransfer coefficient, and external- 
surface temperature at stagnation point. 
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Figure - Diagram of typical chordwise variation of external heat-t rans fe r 

coefficient over propeller blade. 
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(c) Modification 2* Cross-sectional flow ratio, Ap/AjJ = 1.5} Internal 

heat-transfer area ratio, ^ kg - 2.5. 

Figure 12. - Comparison of original with two p rope I I e r- b I ade-sect i on modifi- 
cations that indicate savings in internal heat flow. 
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Surface distance from stagnation point, s 


Figure 13. - Diagram of typical chordwise variation of heat-t ransfer datum 

temperat u re. 
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Fig. 14 





Figure 14. - Comparison of leading-edge surface temperature by two meth- 
ods for station 24^ condition B. 





